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We discuss three mechanical effects initiated by cosmic rays which may limit the 
sensitivity of gravitational wave antennas. Unsolved problems are formulated and 
several recommendations for the antenna designs are presented. 
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q I antennas and antennas based on free masses, e.g 



Since 1969 several groups of researchers have presented results of the analysis and estimates of 
possible contribution of cosmic rays to the noise floor of gravitational wave antennas (both bar 
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The currently achieved sensitivities in 



^ Initial LIGO project (see [Hl^]) and the planned sensitivity in the next stage (Advanced LIGO which 
$h^! is planned to operate within a few years from now) are respectively h ~ 1 0~ 21 and h, ~ 1 0~ 22 for the 
^ amplitude of the perturbation of the metric (at the mean frequency f ~ 100 Hz in the bandwidth 

•i-H . 

^ ■ ~ 100 Hz). Independently of Advanced LIGO several groups of researches which belong to the LSC 
5r (LIGO Scientific Collaboration) are continuing the analysis of new topologies and designs of the 
antennas which will permit them to reach a sensitivity even better than in Advanced LIGO (see e.g. 

)• 



Apart from this activity during the last two decades there is substantial progress in the col 



ection of 



BBS 



data and in the resolution in measuring features of cosmic ray showers (cascades) (see e.g. 
Thus it is reasonable to revise the contribution of cosmic rays to the noise in gravitational wave 
antennas. 

In the present note we limit ourselves to only three possible mechanical "actions" on the rest 
masses (mirrors): 

1. Direct transfer of mechanical momentum from the cascade to the LIGO mirror. 



2. Distortion of the mirror's surface due to the heating by the cascade and subsequent thermal 
expansion — thermoelastic effect. 
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3. Fluctuating component of the Coulomb force between electrically charged mirror and grounded 
metal elements located near the mirror's surface. 

II. DIRECT TRANSFER OF CASCADE MECHANICAL MOMENTUM TO LIGO 

ANTENNA MIRROR 

In Advanced LIGO it is planned to use nearly cylindrical mirrors in the main Fabry-Perot (FP) 
optical resonator (mirror's diameter - 2R ~ 35 cm, height - H ~ 20 cm, total mass of fused silica 
(SiOz) mirror M = 40 kg). The mirrors will be suspended on Si02 fibers in a horizontal vacuum 
tube. The axis of the mirrors will be parallel to the tube. The distance between mirrors will be 
L ~ 4 km. Two mirrors have to respond to the gradient of acceleration generated by a gravitational 
wave propagating in a perpendicular direction to the FP axis. 

TABLE I: The parameters of high energy cascades we used for estimates. S is cascade energy, Jh, J e 
are the fluxes of cascades produced by muons, hadrons and by soft component, consequently, at sea level; 
Ne.max is a number of electrons in the cascade maximum; AS is the energy lost by the cascade in the 20 
cm of S1.O2; N ev is the expected number per year of events with energy losses higher than AS. 
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It is reasonable to expect that a certain part of all cascades (showers) generated by very high 
energy muons of cosmic rays will pass through the mirror along axes which have a small angle with 
the FP resonator axis. 

Thus one may expect that not a very small fraction of a cascade with energy £ will be lost in the 
"travel" through 20 cm of fused silica mirror. This fraction AS will give a momentum AP = AS/c 
to the mirror along its axis (c is speed of light). Correspondingly this mirror position will change 
during time t by a value 

A£t 

Ax ~ — — , t ~ 0.01 s 
Mc 
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The narrow cascades in the direction close to vertical can be produced by so-called unaccompanied 
hadrons j^] and electron-photon (soft) component of cosmic rays High energy (S > 0.5TeV) 

cosmic ray muons are the origin of electromagnetic and nuclear cascades. These particles are born in 
decay of 7t ± - and K 1 * 1 -mesons and charmed particles generated in hadronic showers. The showers are 
initiated by interactions of primaries with the nuclei of air at high altitudes in the atmosphere. The 
energy spectra of all particles (and all cascades generated by them) have power law shape F(£ c ) ~ S~ y , 
y is the power index, y = 2.73 ± 0.05 for [8J; y = 2.5 ± 0.1 for J H and y = 2.2 ± 0.2 for J e |l(|. 
For the estimation of the parameters of cascades it is possible to use a well elaborated model of this 



type of cascades [ll | . 

In Table |U we present numerical estimates for "an appropriate candidate" for the discussed process. 
For several energies S of cascade one can find from the literature the fluxes of cascades produced by 
muons (Jjiund) underground [3,0], hadrons (Jh) Q| and by the soft component (J e ) [10] at sea level 
and calculate the mean number N e>max of electrons in the cascade maximum. The flux of cascades 
produced by muons at sea level (J^) was estimated from experimental data obtained underground 
(J mmd) m taking into account the calculations of energy spectrum and angular distribution of muons 
at sea level The energy lost AS of a cascade was obtained by integration of energy losses of 

electrons in the mirror taking into account the energy spectrum of electrons. The expected number 
N ev of events per year with energy losses equal to or higher than AS was calculated by the formula 
Nev = [J p. x Si + (Jh+ J e ) x x T, (Si, S2 are the areas of the mirror perpendicular and parallel to 
the axis, respectively, T is one year T « 3 x 10 7 s). These estimates were made under the assumption 
that the cascade is coming into the mirror being well developed (more than 6 interactions of the 
fastest electrons). 

The values of AS, presented in Table [Q permit us to estimate the mirror's displacement Ax ~ 
(0.8-7-3) x 1 0~ 18 cm. This numerical value of Ax is less than the amplitude of the sensitivity planned 
in Advanced LIGO: AL ~ HL/2 ~2x 10~ 17 cm. 

Evidently one may expect a much more frequent rate of events: e.g. a cascade with initial energy 
S = 1 GeV will fly through the same mirrors several times per second. But because the value of AS 
in this case will be approximately 3 orders smaller the net effect will be negligible compared to the 
planned sensitivity of Advanced LIGO. 
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III. DISTORTION OF THE MIRROR'S SURFACE DUE TO THE THERMOELASTIC 

EFFECT AND HEATING BY THE CASCADE 

The lost in the mirror's bulk energy of the cascade is likely to be distributed into two parts. The 
first one is the rise of free energy of the solid (creation of new dislocations, of new clusters, etc) and 
the second part produces direct heating of a narrow channel. There is no rigorous analysis in the 
quantum theory of solids which permits us to get a reliable value of the ratio of these two parts. But 
it is very likely that the second one (the heating) will dominate. 

The energy AS is converted into heat over a length H = 20 cm (thickness of mirror). This heat 
will be produced in the trace of the cascade in the volume 7tR 2 H where R c ~ 1 4- 7 cm is the radius 
of the cascade trace. Assuming that the volume ~ R 3 on the mirror's surface can freely expand we 
obtain an estimate of the height AH of the " hill" with footprint ~ R 2 on the surface due to thermal 
expansion: 



R AS 

Here p ~ 2.3 g/ cm 3 is the density, C ~ 7 x 1 6 erg/ g cm 3 K is the heat capacity and <x ~ 5.5 x 1 0~ 7 K _1 
is the thermal expansion coefficient of fused silica. The height AH av averaged over laser beam spot 
with radius r ~ 1 cm is approximately equal to 

R 2 f 2x10- 18 cm if AS = 60 GeV, if R c = 1 cm 
AH av -AHx^~ (3.2) 

r |5.4x10- 17 cm if AS = 230 GeV, if R c = 7 cm 

The displacement AH av of the surface considered above is produced by a cascade developing mainly 
perpendicular to the surface of the mirror. However, there are cascades having traces approximately 
parallel to the surface. Such a "parallel" event produces a greater contribution to the fluctuational 
displacement of the mirror's surface: 

AS R c f 2x10- 17 cm ifA£ = 60GeV 

AH aVjParall ~ x R c a x — ~ i (3.3) 

P CR c 2r r 7.7 x 10- 17 cm if AS = 230 GeV 

However, such events are rare than "perpendicular" ones by a factor of about R 2 /r 2 ~ 0.01 4 0.5 
if we roughly assume the spherical symmetry of showers distribution 1 . 

1 For more accurate consideration we have to take into account the real spartial distribution of muon , hadron 
and soft components of cascades. 
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IV. FLUCTUATING COULOMB FORCE 

In 1995 R. Weiss pinpointed the potential danger from electrical charge accumulated on the 
mirror's surface Q]. Direct measurements of the values of electrical charge density ct on models of 



mirrors were performed independently by several groups 
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demonstrated that the values of ex on models of mirrors (fabricated from fused silica) was from 1 6 
to 1 7 electrons per cm 2 and in several cases even higher. 

Recently V.P. Mitrofanov and his colleagues Q] have measured the values of cr in the same 
vacuum chamber in which record-high quality factors of pendulum and violin modes (Q > 10 s ) 
were demonstrated. These measurements were performed during several months and a slow, long 
lasting drift of a was observed. The high values of surface density of charge obtained mean that the 
electrostatic potential of the mirror may exceed 1 00 V. 

V.P. Mitrofanov and his colleagues [l6| have discovered a monotonic rise of negative charging - 
dcr/dt ~ 10 5 electrons per cm 2 per month. This effect can be qualitatively explained by the model of 
transition effect occurring during transition of cascade particles, soft component and gamma-quanta 
of natural radioactivity through the iron "envelope" to the fused silica mirror (i.e. vacuum chamber). 

In the cosmic rays of low energy there are more electrons than positrons. Compton effect, photo- 
electric effect and the process of ionization are the sources of the excess. In electromagnetic cascades 
the number of particles with low energy is much larger than the number of particles with high 
energy. The number of particles having energy less than 0.05 of the critical energy in the matter 
(£cr,Fe = 20, 7MeV for iron, £ cr>A i = 40MeV for aluminium, £ cr ,sio 2 = 47, 3MeV for SiOi) is equal 
to ~ 20%, and the number of particles having energy less than 0.02 of the critical energy in the 
material is equal to ~ 10% of number of particles in the cascade maximum. The cosmic rays are 
very sensitive even to thin layers of matter. If a cascade developed in the heavy material comes to 
the material consisting of lighter atoms, it brings to the light material more electons than it takes 
away. This is explained by the fact that the number N e of electrons produced in the material is 
proportional to the fraction N e ~ £/£ C r, where £ is the cascade energy. 

As an example we consider a cascade in the maximum of its development coming from iron to 
fused silica. In this case cascade theory gives the formulas for the number of low energy electrons 
produced in iron and in fused silica: 

N e>Fe (h < IMeV) = 0.2 x - ° 2 x J—, (1 MeV = 0.05£ CT , Fe ), (4.1) 

A/ln(£:/E CTiFe ) <W e 

N e , Si02 (E < IMeV) = 0.1 x 03 x -r^— , (1 MeV = 0.021 £ cr , S io 2 ). (4.2) 

^\\l{£/£ cr>si02 ) £cr,Si0 2 
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So we see that the number of low energy electrons produced in iron and coming to the mirror is 
about 3 times larger than the number produced in fused silica and outgoing from it: 



N 



e,Fe 



0.4 S CT 



SiQ 2 



Ne,Si0 2 0.3£ cr> Fe 

This electron excess will stay near the surface of the mirror and will give an additional charge 
to it. The estimates of the number of electrons with energy less than 1 MeV coming from iron to 
fused silica are given in the Ta ble ITT1 This effect can qualitatively explain the monotonic rise of 
negative charge observed in jl^l . 
be performed. 



17| . However, for quantitive explanation a detailed analysis has to 



TABLE II: The mean number N of electrons with energy less than 1 MeV coming from iron to fused silica. 
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The initial design of the "entourage" of the suspended mirror includes several parts which are 
planned to be made of metal. These parts include a "cradle" situated under mirror (this " cradle" 
has to catch the mirror if one or all fibers break). Other parts are called "stoppers" which have to 
limit large horizontal swings of the mirror in case of an earthquake. These parts have to be grounded. 
Thus due to electrical charging of the mirror it is very reasonable to expect that a d.c. Coulomb 
force may act on the mirror. If a grounded metal part has a flat surface S that is close to a part of 
the mirror's surface, then 

F dc ~ InSo 1 ~ 1 .5 x 1 CT 2 dyn, if S = 1 2 cm 2 

This numerical estimate of the d.c. force has to be taken into account in the design of feedback 
actuators which have to maintain the distance between antenna mirrors with accuracy better that 
\/T ~ 10~ 9 cm [T is the finesse, A is the optical wave length). More important is another effect: 
The a.c. component of the Coulomb force may mimic the force F grav that antenna has to register: 



F ~ 

1 grav — 



HLMuj 2 v 7 
— ^ ~ 3 x 1(T 7 dyn 



(4.3) 



A comparison of the values presented above indicates that a relative fluctuation Acr/cr ~ 10~ 5 of 
surface charge density will inevitably produce a " step" of F ac approximately equal to the amplitude 
Fgrav which is the goal of Advanced LIGO. Note that N£<iMeV electrons outgoing from iron to fused 
silica presented in Table ITT1 comes to the square about TtR 2 ~ 100 cm (R c is the radius of cascade) 
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. Then one can estimate the relative fluctuations of charge density Act/ ct caused by a single cascade 
with energy £ = 2 TeV: 

— ~ 1Cr 6 ^2 x 1(T 5 

CT 

We see that this fluctuation is strong enough to produce an a.c. component of Coulomb force larger 
than F grav if a relatively large surface of grounded metal plate will be located near the mirror. 

V. CONCLUSION 

It is evident that the first two effects, being not very strong ones for Advanced LIGO, may be 
relatively easily vetoed by requiring coincidence between detectors (if at least two antennas are 
operating). On the other hand, the veto can not be considered as an absolute "remedy" for low 
values of the signal to noise ratio. It is worth noting that in the next stage after Advanced LIGO 
these two effects will make serious contributions to the level of the noise floor. 

This conclusion may not be automatically extended to the third effect. First of all because the 
negative charging may be high when the mirrors "spend" a long time in the vacuum (an year or 
longer) without scheduled removal of the accumulated electrical charge. The second "reason" is the 
design of all metal parts of the mirror's "entourage". One recommendation for this design is evident: 
it is necessary to use small square areas of all metal element that are near the mirror's surface and 
to place these elements as far away as possible. 

There are two evident recommendations for the consequent measurement and analysis: 

1. To measure bursts of electrons which appears on the mirror's surface with resolution better 
than 10 2 e/cm 2 and time shorter than 10~ 2 sec. 

2. To analyze the possibility to cover the mirror's surface over the coating with a transparent few 
nanometers thick layer with substantial conductivity to reduce the d.c. component of electrical 
charge. 

In all three effects considered above the mechanical action on the mirror produces a step-like 
displacement (either of the mirror's center of mass or of its surface) . This type of response is similar 
to the one predicted for the shape of gravitational wave bursts created in the process of supernova 



explosion predicted by V.Imshennik 18|. The predicted rate of these bursts is approximately two 



orders higher than the rate of neutron star merger events. 
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